Perturbation Bounds for Isotropic Invariant Subspaces of Skew-Hamiltonian Matrices
نویسنده
چکیده
Abstract. We investigate the behavior of isotropic invariant subspaces of skew-Hamiltonian matrices under structured perturbations. It is shown that finding a nearby subspace is equivalent to solving a certain quadratic matrix equation. This connection is used to derive meaningful error bounds and condition numbers that can be used to judge the quality of invariant subspaces computed by strongly backward stable eigensolvers.
منابع مشابه
Perturbation bounds for $g$-inverses with respect to the unitarily invariant norm
Let complex matrices $A$ and $B$ have the same sizes. Using the singular value decomposition, we characterize the $g$-inverse $B^{(1)}$ of $B$ such that the distance between a given $g$-inverse of $A$ and the set of all $g$-inverses of the matrix $B$ reaches minimum under the unitarily invariant norm. With this result, we derive additive and multiplicative perturbation bounds of the nearest per...
متن کاملRelative perturbation bounds for the eigenvalues of diagonalizable and singular matrices – Application of perturbation theory for simple invariant subspaces
Perturbation bounds for the relative error in the eigenvalues of diagonalizable and singular matrices are derived by using perturbation theory for simple invariant subspaces of a matrix and the group inverse of a matrix. These upper bounds are supplements to the related perturbation bounds for the eigenvalues of diagonalizable and nonsingular matrices. © 2006 Elsevier Inc. All rights reserved. ...
متن کاملPerturbation Theory for Simultaneous Bases of Singular Subspaces
New perturbation theorems are proved for simultaneous bases of singular subspaces of real matrices. These results improve the absolute bounds previously obtained in [6] for general (complex) matrices. Unlike previous results, which are valid only for the Frobenius norm, the new bounds, as well as those in [6] for complex matrices, are extended to any unitarily invariant matrix norm. The bounds ...
متن کاملNumerical Computation of Deflating Subspaces of Skew-Hamiltonian/Hamiltonian Pencils
We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear-quadratic optimal control problems, H∞ optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requires computing invariant and deflating subspaces of matrices and matrix pencils with Hamiltonian and/or ske...
متن کاملNumerical Solution of Quadratic Eigenvalue Problems with Structure-Preserving Methods
Numerical methods for the solution of large scale structured quadratic eigenvalue problems are discussed. We describe a new extraction procedure for the computation of eigenvectors and invariant subspaces of skew-Hamiltonian/Hamiltonian pencils using the recently proposed skew-Hamiltonian isotropic implicitly restarted Arnoldi method (SHIRA). As an application we discuss damped gyroscopic syste...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM J. Matrix Analysis Applications
دوره 26 شماره
صفحات -
تاریخ انتشار 2005